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Abstract. - Considering a two-component Bose-Einstein condensate in a double- well potential, 
a method to generate a Bell state consisting of two spatially separated condensates is suggested. 
For repulsive interactions, the required tunnelling control is achieved numerically by varying the 
amplitude of a sinusoidal potential difference between the wells. Both numerical and analytical 
calculations reveal the emergence of a highly entangled mesoscopic state. 
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In a recent ground-breaking experiment self-trapping 
was observed for small Bose-Einstein condensates (BECs) 
m a double- well potential [1]. A major experimental goal 
in this system is the generation of mesoscopic entangled 
number states like "Schrodinger-cat" states [1]. Unfortu- 
nately, in a double well the ground state of a condensate 
with repulsive interaction is not the desired "Schrodinger- 
cat" state (e.g. ref. [2]). On the other hand using attrac- 
tive condensates leads to the problem of instability (e.g. 
ref. [3]). However, a two-component repulsive condensate 
in a double well is investigated in this letter, because it 
(Offers the possibility to increase the entanglement of the 
ground state drastically. 

] To control the entanglement generation, we employ a 
method recently suggested to induce the phase transition 
between a superfluid and a Mott insulator [4] by apply- 
ing time-periodic potential differences [5,6]. While on the 
single particle level the effects predicted for periodically 
driven systems [7, 8] have not yet been realised experi- 
mentally, the coherent control by changing the amplitude 
of periodic force fields might be easier to achieve for BECs 
rather than for single particles. 

Suggestions to produce mesoscopic entangled states, in- 
volve coherent scattering of far detuned light fields [9]; 
several groups suggest to produce mesoscopic entangle- 
ment with condensates by manipulating the interaction 
between the particles, or by controlling the dynamics of 
the system [10-15]. While the experimental realisation of 
such mesoscopic entangled states is still a challenge of cur- 
rent fundamental research, possible applications include 
quantum computing or quantum information processing. 

In this letter we study a double-well potential similar 
to the experiment [1], but with a two-component conden- 



sate (for various aspects investigated for two-component 
BECs see e.g. refs. [16-19]). To motivate the mechanism 
by which we control the entanglement generation, we start 
with the two-particle case, because it exhibits all crucial 
features, before we extend both the numerical and analyt- 
ical calculations to N particles. 

The model. — We study a symmetric double-well po- 
tential filled with Bose particles at very low temperature. 
Adopting the common two-mode approximation, assum- 
ing only on-site interactions and denoting the tunnelling 
splitting between the two lowest single particle energy 
states with hQ^ we use a well established model describing 
the dynamics of a one-species BEC in double-well poten- 
tial [20]. This model gives a good qualitative description 
of recent experiments [1]. Extending the model to a bi- 
nary condensate consisting of N particles of species A and 
N particles of species B leads to the Hamiltonian 
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where a-^^ and 6-^^ are the annihilation (creation) oper- 
ators for a boson of species A or B in the ith well (i = 
1,2) and h/jjf{t) specifies an externally applied potential 
difference between the two wells. The relation between 
the interaction parameters hca, i^b and k^ab and the ex- 
perimental situation can be found e.g. in ref. [2]. We use 
the Fock-basis N — i)A\j, N — j)^ = where i/j are 
the numbers of A/B particles occupying the first well and 
{N — i)/{N — j) the number of A/B particles occupying 
the second well. 
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In these notations, the Bell state reads: 

VBell = -^(|iV,0) + |0,iV)) 

= -^(|iV,0)A|0,iV)B + |0,iV)A|7V,0)B) . (2) 

In this state, the two condensates are maximahy entangled 
with each other in the sense that each condensate is in a 
superposition of being in both wehs and if one condensate 
is measured in one weU, the other condensate wiU be in 
the other well. 

Similar spatially entangled states are known from the 
EPR-paradox and are of special interest in the field of 
quantum cryptography, where especially entangled pho- 
tons are used [21,22]. Looking at the AT-particle dynamics, 
we have to take (A/" + 1)^ states into account. If we exem- 
plarily study the case of only two distinguishable particles 
A and B without driving the well (/i = 0), the system can 
be solved analytically. While we assume short range inter- 
action between the two particles typical for cold atoms, the 
resulting entanglement is similar to the case of two elec- 
trons in a double quantum dot [23]. Here, we use the basis 
{|1, 1), |1, 0), |0, 1), |0, 0)} to obtain the system's Hamilto- 
nian as a 4 x 4 matrix 



Bell state by just generating a binary-BEC in a very deep 
double- well potential; the system could be in any super- 
position of those two states. If, however, one starts in the 
ground state of a not so deep well, the Hamiltonian does 
not allow transitions between symmetric and antisymmet- 
ric states (cf. ref. [24]). Thus, thermal excitations to the 
(antisymmetric) first excited state can be discarded. 

The double-well potential is now modulated periodi- 
cally in time by the force / indicated in eq. (1). As 
the new Hamiltonian is periodic in time, Floquet theory 
can be applied. As shown in refs. [25,26], for sufficiently 
high frequencies {uo/^ ^ 1) the system behaves like an 
undriven system with a rescaled effective tunnelling fre- 
quency l^eff = ^ • Jo (2/i/c(;) or equivalently with the effec- 
tive interaction parameter 
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where u = a^ab/^ is the interaction parameter. When 
repulsive interaction is assumed ('ix > 0), the ground state 
and first excited state of the system are given by 



where Jo denotes the ordinary Bessel function of order 
zero. The system is coherently controlled by driving it 
with linearly in time increasing parameter 2/i/cj (in this 
case 2ii/uo 2.405 • r/rmax, where r = Q • t is a dimen- 
sionless time variable). At the time Tmax the parameter 
reaches the first zero of Jo; the effective interaction be- 
comes very large. Similar effects could be achieved with- 
out driving by adiabatically increasing the depth of the 
wells. Compared with this coherent control mechanism, 
the feed-back loops of optimal control theory [27] would al- 
low to further optimise entanglement generation for given 
experimental parameters. 

To quantify the entanglement of a state |^) we use the 
fidelity, i.e. the probability of the system to be in the Bell 
state given by the square of the state's overlap with V^Beii- 
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ing energies ^o = hVt{u/2 — 

the case of a very large interaction parameter the ground 
state approaches the Bell state (eq. (2) with TV = 1), 
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and becomes a maximally entangled state. On 
the other hand with no interaction {u = 0) the 
ground state is given by 1/2 .(1111)= 1/2 • 
(|1,1) + |0, 1) + |1,0) + |0,0)), which can be written as 
a product state yo{u = 0) = 1/2 (|1, 0)a + |0, 1)a) • 
(|1, 0)b + |0, 1)b) and therefore is not entangled at all. 
Furthermore, for very large interaction parameter u the 
ground state is quasi- degenerate. This is also true for the 
TV-particle system with the states (|iV,0) =b \{),N))/^. 
Because of this quasi-degeneracy one cannot generate a 



Because our goal are fidelities close to one, more sophisti- 
cated entanglement measures (see e.g. refs. [22,28]), which 
identify much less entangled states, are not necessary. As 
can be seen in fig. 1 entanglement is maximized by the 
procedure described above when the amplitude is changed 
adiabatically (for a value of say Tmax = 2 the increase 
would be too fast and the maximal entanglement is not 
reached). On the other hand the driving frequency must 
+ 1 1,0)) , (5) ^ot be too low. For uj/Q. = 100 everything works excel- 
lently as expected; even ioi uo/Vt = 10 we get very good 
results, but with uo/VL on the order of one a monotonic 
increase of Peeii is not observed. 



Analytical calculation for the Bell state. — In the 

spirit of refs. [2,10] we use the mean-field equations, in or- 
der to extend the two-particle dynamics to many Bosons, 
and construct a symmetrised superposition of two mean- 
field states. Within the mean-field approximation the 
system is reduced to four amplitudes describing the dy- 
namics of condensate A (B), such that |aip = 1 — |a2p 



1 — 1^2! ) is the expected fraction of condensate A 
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Fig. 1: Fidelity (eq. (7), probability of the system to be in the 
Bell state) for the system with N = 1 particle of each species. 
After initialising it in the ground state the driving parameter 
2ijl/uj was linearly increased in the time r until reaching the 
first zero of Jo(2/i/a;) with (2/i/a;)max — 2.405; the interaction 
parameter was u = kab/^ = 0.5. The driving frequency used 
is Lu/Q = 10 (upper curve) and lu/Q = 1.2 (lower curve). The 
increase with the higher frequency is in almost perfect agree- 
ment with the ideal curve indicated by the dashed line. For the 
lower frequency the increase does not work and the Bell state 
is not reached. 



(B) in the first well. The equations of motion for two con- 
densates with equal number of particles (TVa = = 
are given by 
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where the dot means the derivative with respect to r = ^l-t 
and ka, a^b and kab are the interaction parameters given 
in eq. (1). In the following the intra-condensate interac- 
tions are equalised and the abbreviations a = 2^ and 
/3 = 2 ^^^^ are used. Experimentally, a similar situation 
could be achieved by choosing the numbers of particles in 
the condensates, such that Nak.a = A^b^b- The inter- 
condensate interaction parameter f3 can be modulated by 
a Feshbach-resonance [29,30]. 

To find stationary state solutions we apply the ansatz 

(^1,2 = '01,2 exp(— ii/t) and bi^2 = 0i,2 exp(— i/it) . (12) 

One finds two kinds of solutions, a balanced one, where 
the condensates are equally distributed on the two wells 



m = v^i = 0' 



I) and an unbalanced solution. 



where condensate A/B is mainly in the first/second well. 
For {2f3 — a) > the amplitudes in this case are 
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A second unbalanced solution is obtained by interchang- 
ing the indices 1 and 2. If (2/^ — a) is negative, ipi and 
^p2 in equation (13) have opposite signs. As the resulting 
state is not the ground state in the limit of large interac- 
tions necessary for entanglement generation, this solution 
can be discarded here. The mean- field Hamiltonian corre- 
sponding to eqs. (8)- (11) reads: 



H 



[-{aia2 + ala2) + ce(|ai|^ + l^^2|^)] 
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The energies for the balanced and unbalanced states per 
particle therefore are given by 
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When is the unbalanced state the ground state? It is, if 
^ubai < ^bai, leading to 



a- p-\-2- 
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When again by driving the double-well the effective in- 
teractions are strongly increased (6) only the difference of 
the intra- and inter- condensate interactions a and P de- 
termine whether the balanced or the unbalanced state is 
the ground state. 

To compare the nonlinear mean-field dynamics to the 
A^-particle dynamics, we need the TV-particle counterparts 
of the mean-field states. By using the atomic coherent 
states [2] 



^ \nA 

nA=0 ^ ^ 



cos^^(7^/2) sin^-^^(7^/2) 

(17) 



for condensate A (and B analogously) with cos(i?/2) = 
-01, sin(7^/2) = i\)2 and phase 99 = we can construct a 
symmetrised superposition of two unbalanced states 



I*) = M 0)a|i?, 0)b + 0)a|i?, 0)b) 



(18) 



where the overlined term is obtained by replacing 
|n,A^-n)A/B by \N - n,n)A/B- The probability Peeii 
is then given by 

PBeii = 2ArMV'r+^l'^)' • (19) 
For large TV the above expression is approximated by 

„pf-S^V (20) 
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Fig. 2: Fidelity (7) for the system with = 50 particles of each 
species (cf. fig. 1). After initialising it in the ground state the 
driving parameter 2/i/a; was linearly increased until reaching 
the first zero of Jo(2/i/a;) ; the used driving frequency was 
uj /Q. — 40. The upper plot shows that generating the Bell state 
is possible on a short time scale with a fidelity higher than 
96%; the interaction parameters were ka/^ — n^/^ = 0.1 
and a^ab/^ = 0.115. Inset: If the interaction parameters are 
chosen such that the initial state is not bimodal (see fig. 3) 
with ka/^ = kb/^ — 0.1 and /^ab/^^ = 0.11, it takes longer 
to adiabatically reach the Bell state. The fidelity is again very 
high (> 97.5%). Lower figure: although the unbalanced state 
is already quite entangled when it becomes the ground state 
for r ^ 80 (cf. eq. (16)), the analytic curve indicated by the 
dashed line (cf. eq. (20)) gives a good qualitative description of 
this process. The used interaction parameters here are ha/^ — 
= 0.0016 and ha^/Q. = 0.0133; the fidelity is > 97.5%. 



TV-particle dynamics and entanglement genera- 
tion. — To demonstrate entanglement generation numer- 
ically, we use 100 particles (50 of each kind), which is close 
to the experimentally reachable domain [31]. In experi- 
ments the decision for small condensates makes sense, be- 
cause decoherence times induced by particle loss are long. 
Again we prepare the system in the ground state and drive 
it with linearly in time r increasing driving parameter 
2/j./ijj up to the first zero of Jq. Figure 2 shows the emer- 
gence of the highly entangled state under the influence of a 
periodic force-field. The parameters are chosen depending 
on the desired behaviour: generation of entanglement on 
very short time scales (fig. 2, top), good agreement with 
the analytical calculation (eq. (20)) or generation of entan- 



Fig. 3: Probabilities of the system to be in the Fock state |i, j), 
where i and j are the Fock indices. The system was initialised 
in the ground state, which is in this case (same parameters 
as in the inset of fig. 2) a quite balanced state (upper figure: 
system at r — 0), and then driven with linearly increasing 
driving parameter 2fi/uj. With increasing time the system's 
state becomes more and more unbalanced (lower figure: system 
at r = 250) until it reaches the Bell state at Tmax = 500 (cf. 
inset of fig. 2). 



glement starting with a balanced initial state (see fig. 3). 
Although the mean-field equations are valid in the limit 
N ^ oc and k ^ such that Nk is constant, the predic- 
tions obtained by our calculations qualitatively agree with 
the numerical results as long as (2/5 — a) is not too high in 
the regime, where the balanced state is the ground state. 

In fig. 3 we show the occupation of states by plotting 
the Fock occupation numbers i and j and the correspond- 
ing probabilities The probabilities move to the 
outer (more unbalanced) states while increasing the pa- 
rameter 2fi/u; and at the end only the two states |iV, 0) 
and |0, A/') are nearly equally occupied. So a quite bal- 
anced state is changed into a totally unbalanced and max- 
imally entangled state. 

Within the validity of the two-mode approximation (cf. 
ref. [15]), the fidelity (7) can be raised by using higher 
driving frequencies uj/Q. As the entanglement generation 
involves adiabatic following of Floquet states (see ref. [5] 
and references therein), the time-scales on which the am- 
plitude will have to be changed in future experiments can 
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Fig. 4: Part of the quasienergy spectrum for the system (1) 
with N = 10 particles of each kind and driving frequency 
uj/n 100, k.a/^ = f^B/ft = 0.1 and k.b/^ 0.125. 
At 2/i/a; 2.4 the ground-state, the lowest line, is quasi- 
degenerate (as many other states, cf. ref. [5]). At this point 
Jo(2/i/(x;) is nearly zero and therefore tunnelling is strongly 
suppressed. The coherent control used in this letter involves 
adiabatic following of the ground state. While transitions to 
the state which follows the first excited state are suppressed 
for symmetry reasons (cf. ref. [24]), the energy difference to the 
higher states gives a measure for the timescale on which the 
amplitude can be changed. (For N = 50 most of the spectrum 
would become dense. However, for high driving frequencies the 
lowest two quasienergies are still well separated from the other 
quasienergies - otherwise the high fidelities of fig. 2 would not 
have been reached). 



be deduced from calculating the quasienergies (fig. 4). For 
higher particle numbers than chosen in fig. 4 and not too 
high driving frequencies, avoided crossings occur which are 
not important for entanglement generation as long as the 
driving frequency still is in the high frequency limit (fig. 2; 
see also ref. [5]). To find the best parameters for a given 
experiment with respect to low decoherence and large en- 
tanglement and to optimise the protocol /i(r), i.e., the 
way the driving amplitude is changed in time, one might 
employ the methods of optimal control theory [27]. 

Conclusion. — We demonstrated that the generation 
of a highly entangled state with a binary condensate in a 
double- well potential is possible. By driving the double- 
well with an increasing amplitude the ground state is 
turned into a Bell state provided that the inter-condensate 
interaction is first slightly increased via a Feshbach- 
resonance. The fidelities of our results within the model 
of eq. (1) are even on short time scales above 96% which 
is considerably higher than the fidelity obtained in similar 
models for single-species condensates [13-15]. 
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